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The main contribution of this study is to open a discussion 
regarding the structural optimization associated with the cost 
efficiency and structural reliability sufficiency consideration. To do 
so, several various optimization approaches are investigated to 
deliberate both cost and reliability concerns. Particularly, particle 
swarm optimization is highlighted as a_ reliable optimization 
approach. Accordingly, an illustrative example is rendered to 
compare the feasibility of the considered optimization approaches. 
The feasibility of the investigated approaches is evaluated using the 
cost and reliability analysis. For the considered example, it was 
observed that the PSO optimization algorithm has _ multiple 
advantages such as_ easy realization, fast convergence, and 
promising performance in nonlinear performance optimization. The 
PSO optimization algorithm can be successfully applied in various 
fields of civil engineering. This popularity is due to the 
understandable performance of the PSO as well as its simplicity. In 
this paper, first, the literature on the subject has been described by 
two-dimensional truss analysis using the finite element method and 
optimized using the PSO particle swarm algorithm. A comparison 
of the results with this reference indicates the accuracy of this 
particle swarm algorithm in truss optimization. Indeed, this study 
ignites two main insights in structural optimizations assessment. 
The first illustration is related to how to establish a framework for 
structural system reliability analysis associated with the different 
degrees of indeterminacies. And the second illustration is related to 
making a decision problem concerning the structural optimization 
while both cost and reliability metric are two main parameters for 
the construction point of the view. 


How to cite this article: Rezaee Manesh M, Ghasemi SH, Rezaee Manesh M. Dual target optimization of two-dimensional truss 


using cost efficiency and __ structural 


reliability sufficiency. J Soft Comput Civ Eng 2020;4(4):98-111. 


https://doi.org/10.22115/SCCE.2020.244833.1252. 


2588-2872/ © The Authors. Published by Pouyan Press. 


This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/). 


M. Rezaee Manesh et al./ Journal of Soft Computing in Civil Engineering 4-4 (2020) 98-111 99 


1. Introduction 


Swarm intelligence (SI) is the collective behavior of decentralized systems consisting of many 
individuals who coordinate their activities using self-organization [1]. The SI systems are 
inspired by biological systems that are usually composed of a population of simple agents 
interacting with each other and globally with their environment. Ant colonies [2], bird flocks 
[3,4], honey bee colonies [5,6], and fish schools [7,8] are examples of natural SI systems. Based 
on the systems, several optimization algorithms such as ant colony, artificial honey bee colony, 
and particle swarm optimization (PSO) algorithms have been developed. PSO is a computational 
method based on a powerful population that is capable of solving the problem using the 
population of candidate solutions. To improve the candidate solution, PSO frequently moves the 
candidates using simple mathematical equations in a search space [9]. Since the introduction by 
Kennedy and Eberhart (1995) [10] and Eberhart and Kennedy (1995) [11], the PSO algorithm 
has been successfully used in various fields of civil engineering. 


The optimization methods and algorithms are divided into two categories: exact algorithms and 
approximate algorithms. Exact algorithms can accurately find the optimal solution, but they do 
not work in hard optimization problems, and the solving time in these problems increases 
exponentially. Approximate algorithms can find good (quasi-optimal) solutions to difficult 
optimization problems in the short solving time. Approximate algorithms are also divided into 
two categories: heuristic and meta-heuristic algorithms. The two major problems with heuristic 
algorithms are (1) falling in local optima, and (2) inability to be applied to different problems. 
The meta-heuristic algorithms have been proposed to solve the problems of heuristic algorithms. 
The meta-heuristic algorithms are one of the approximate optimization algorithms that have the 
mechanism to exit local optima and can be used in a wide range of problems. In recent decades, 
various categories of this type of algorithm have been developed. Over the past few years, PSO 
has been widely used in various aspects of geotechnical engineering such as slope stability 
analysis, pile and foundation engineering, rock and soil mechanics, tunneling, and underground 
space design [12]. 


2. Particle swarm optimization algorithm 


In many optimization problems, especially the big ones, choosing the best solution through the 
global search, though not being impossible, is very difficult. The goal of the optimization 
problem is to reduce the search time. Heuristic methods are good solutions for finding the 
optimal solution, but they do not guarantee to find the optimal solution. Today, however, as the 
problems become bigger, the popularity of heuristic and meta-heuristic methods has considerably 
grown. 


Particle swarm optimization (PSO) is one of the meta-heuristic optimization techniques that 
work based on the population. The original idea for the method was first proposed in 1995 by 
Kennedy and Eberhart [10] which was inspired by the collective behavior of fish and birds for 
finding the food. A group of birds and fish look for food in random space, and there is only one 
piece of food, and none of the birds knows the location of the food and only knows the distance 
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to the food. One of the best strategies is to follow a bird that is closest to the food. This view is 
the main strategy of the algorithm. 


In this method, each bird has a possible solution in the problem space called a particle. Each 
particle has a fitness value that is calculated by the problem. The particle that is closest to the 
solution has the highest fitness. This algorithm has a continuous nature and has been proven to 
be efficient in numerous applications [13]. 


This algorithm is one of the population-based parallel search algorithms that start with a group of 
random solutions (particles) and then continues to search for finding the optimal solution in the 
problem space by updating the position of the particles. Each particle is denoted in the 
multidimensional mode (depending on the type of problem) by two vectors Xid and Vid, which 
represent the position and velocity in dimension d of particle i, respectively. At each stage of the 
population movement, the position of each particle is updated by two best values. The first value 
is the best experience ever gained by a particle and is represented by p_best. The second value is 
the best experience gained among all the particles and is shown by g best. In each iteration, after 
finding these two values, the algorithm updates the new velocity and position of the particle- 
based on Equations (1) and (2). 


vig(t + 1) = w. vig(t) + cy.rand, (p_bestia (t) — xiq(t)) + C2.rand, (9 _bestig(t) - xia(0)) (1) 


Xiq(t + 1) = xiq(t) + via(t + 1) (2) 


In Equation (1), w is the inertia weight factor, which is usually in the range [0-1]. C1 and C2 are 
the learning or acceleration factors, which are selected in the range [0-2], which, in most cases, 
are set to 1.49 or 2. rand1 and rand2 are also two random numbers in the range [0-1]. Also, the 
final value of each particle's velocity to prevent algorithm divergence is limited to one range. 


Via © |—Vmax -Vmax] is the termination condition of the convergence algorithm to a certain 
extent or to stop after a certain number of iterations. Equation (2) also updates the current 
position vector of the particle according to its new velocity. 


The right side of Equation (1) consists of three parts; the first part is a factor of the current 
velocity of the particle, the second part is used for changing the velocity and rotation of the 
particle towards the best personal experience, and the third part is used for changing the velocity 
and rotation of the particle towards the best group experience. 


The optimization of the particle swarm without the first part of Equation (1) will be a process in 
which the search space is gradually reduced and the local search is formed around the best 
particle. In contrast, the first part of Equation (1) will cause the particles to move in their normal 
direction to reach the wall of the search area and perform some sort of global search. By 
combining the two factors, it was attempted to strike a balance between local and global search. 
To better strike such balance, w was first proposed in [14], which specifies the movement factor 
in the global search. The parameters cl and c2 also specify the movement factor in local search. 
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Fig. 1. PSO flowchart [15]. 
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Fig. 2. Movement of a particle in PSO. 
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Current position 


Fig. 3. Schematic structure of particle in PSO [16]. 


2.1. PSO algorithm in unconstrained problems 


The unconstrained PSO algorithm can be expressed as the figure below. 


Initialize particles with random position and velocity vectors 


For each particle’s position “P” evaluate fitness 
Compare particle’s position (P) with fitness(pbest) 


Loop until all particles 


suoneioy xeuw [UN doo7T 


Set best of pbest as gbest 
update particle’s velocity and position 
Slop giving gbest as the optimal solution 


Fig. 4. Flowchart of particle swarm optimization algorithm in unconstrained problems [13]. 


2.2. PSO algorithm in constrained problems [13] 


In the constrained problems, the particles should satisfy the constraint, which is the case with the 
truss problems. 
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Particle swarm optimization algorithm in constrained problems: 


For each particle { 
Do { 
Initialize particle 
} While particle is in the feasible apace (1.e. it satisfies all the constraints) 
} 
Do { 
For each particle { 
Calculate fitness value 
If the fitness value is better than the best fitness value (pBest) 
in history AND the particle is in the feasible space, set current 
value as the new pBest 
} 
Choose the particle with the best fitness value of all the particles as 
the gBest 
For each particle { 
Calculate particle velocity according to equation (a) 
Update particle position according to equation (b) 
} 


} While maximum iterations or minimum error criteria is not attained 


3. Truss analysis using finite element method 


Trusses include two-force members, and the external forces are only applied to the nodes. 


wo <4 
uD < 


SF 
F 


Fig. 5. Two-dimensional truss. 
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3.1. Local and global coordinates 


Figure 6 shows a truss member and the equations for transforming two local and global 
coordinates are derived: 


Deformed 
element 


x 
Fig. 6. A member of truss and equations for transforming two local and global coordinates. 

Local coordinate system 

q' = [a1 -9]" (3) 

Global coordinate system 

Q = [Q1.Q2- Qn ]” (4) 


The relationship between q and q' is as follows: 


qi = 91 cosO + qzsind (5) 
do = q3.cosO + qusin@ (6) 
q =Lq (7) 


Transformation matrix 


_fl m0O0 
a k 0 lm (8) 
Calculation of m and | 
l= cos@ = 2 (9) 


e 


Mm = cosp = sind = om (10) 


e 
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le = VJ (%2 — %1)? + G2 — 91)? (11) 
Element stiffness matrix 


Element stiffness matrix in the local coordinate system 


ki = ete] (12) 


le 
Element stiffness matrix in the global coordinate system 
The formula to calculate the strain energy in the local coordinates can be expressed as follows: 
Ue =5q'7k'q (13) 
q’ =Lq > U. =5 aq" [IT K'L]q (14) 
The formula for calculating the strain energy in the global coordinates is as follows: 


Ue =-q'kq (15) 


2 


Element stiffness matrix in the global coordinate system 


Fea Od "a B (16) 
_fl m OO 
L=[) 0 Im!” os) 
2 Im —-l -lm 
EeAe| lm m2 —Ilm —m? 
k = 22 18 
l | -—l4 —Im F- lm ve) 
—lm — m2? lm wm? 


The stiffness matrices of the elements combine to form the stiffness matrix of the structure. 


Stress calculations 


o = Eye = E.Bq' (19) 

=F [1 {Ul =o (20) 
e 2 e 

q =lq>o==[-1 1Lq (21) 

_fl m0O0O 

=[o 0 lm o 

o=“[-l —m 1 m] (23) 


le 
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4. Numerical example of single-objective truss optimization 
4.1. Particle swarm optimization 


In the example of the optimization of a conventional truss [17], which was analyzed by various 
methods so far, it is optimized by particle swarm optimization (PSO) algorithm method and the 
final result is compared with other methods. For this purpose, the total number of particles is 
considered to be 60, and the constant values of w, ci, and cz are selected as 0.5, 1.5, and 1.5, 
respectively. The algorithm code is written in the MATLAB software and the structure is 
analyzed using the stiffness matrix (finite element) method. 


Wy w/ 
p p 
Fig. 7. 10-bar truss from [17]. 
Table 1 
Parameters related to Fig. 7. 
Parameter Value Unit 
E 104 ksi 
p 0.1 lbm/in? 
Callow 25000 psi 
Cingy 2 in 
L 360 in 
P 100 kip 


Consider the 10-bar truss in Figure 1 which was studied to compare and verify the algorithm. 
The purpose of the problem is to minimize the weight of the structure f(x): 


POH > 1 OAz be) (24) 


x 1s the solution to the problem, Aj is the cross-sectional area of bar i, Lj is the length of bar i, and 
p is the weight density: 


k 
p=0.10 !b/, , (2770 “9/ 3) (25) 
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The elastic modulus FE =1%x10* ksi (6.89 x 10* Mpa) and downward external force 
100 kips (445.374 KN) are applied to nodes 2 and 4. 


The truss constraints are: 


6; <q, i= 1:10 (26) 
Uj < Uq (27) 
where oi is the stress in bar i and 6, is the maximum allowable stress for all bars, u; is the 


displacement of each node (horizontal and vertical), and u, is the maximum allowable 
displacement for all nodes. 


Uq = 2in (50.8 mm) (28) 
Oat = £25 Ksi (172.25 Mpa) (29) 


This problem is optimized using the PSO algorithm, and the results are summarized in Table 2; 
To compare with the results of Table 3, the problem is optimized twice (the cross-sectional area 
of the bars is considered in two different ranges): 


Table 2 
Results from the optimization of 10-bar truss size using PSO algorithm method, first run. 
0.1 < A; (in?) < 33.5, Iteration: 50 
Method Weight (Ib) | Al | A2 | A3 A4 | A5 | A6 | A7 A8 AD | Al0 


PSO 
(obtained Results) 5187.4 31.18 | 0.1 | 26.16 | 16.80 | 0.1 | 0.1 | 7.28 | 20.85 | 21.03 | 0.1 


Table 3 
Results from the optimization of 10-bar truss size using PSO algorithm method, second run. 
1.62 < A; (in?) < 33.5, Iteration: 50 
Method Weight (Ib) | Al | A2 A3 A4 A5 | A6 | AZT A8 AY | Al0 


PSO 
(obtained Results) 5636.2 32.2 | 1.62 | 25.89 | 17.08 | 1.62 | 1.62 | 8.37 | 22.56 | 21.52 | 1.62 


Table 4 
Results from the optimization of 10-bar truss size from [17]. 
Method Weight Al A2 A3 A4 A5 A6 A7 A8 A9 Al10 
(Ib) | (in?) | (in?) | (in?) | (in?) | (in?) | (in?) | (in?) | (in?) | (in?) | (in?) 
OPTDYN | 5472.00 | 25.70 | 0.10 | 25.11 | 19.39 | 0.10 | 0.10 | 15.40 | 20.32 | 20.74 | 1.14 
CONMIN | 5563.00 | 25.20 | 1.89 | 24.87 | 15.83 | 0.10 | 1.75 | 16.67 | 19.73 | 20.98 | 2.51 
GENETIC | 5620.08 | 33.50 | 1.62 | 22.00 | 15.50 | 1.62 | 1.80 | 14.20 | 19.90 | 19.90 | 2.62 
Rajeev 5613.84 | 33.50 | 1.62 | 22.00 | 15.50 | 1.62 | 1.62 | 14.20 | 19.90 | 19.90 | 2.62 
M-3 5719.00 | 25.84 | 3.07 | 26.42 | 12.77 | 0.10 | 3.43 | 19.34 | 19.17 | 18.76 | 4.42 
M-5 5725.00 | 25.83 | 2.88 | 26.45 | 12.75 | 0.10 | 3.77 | 19.37 | 19.18 | 18.77 | 4.38 
GRP-UI_ | 5727.00 | 24.78 | 4.17 | 24.78 | 14.45 | 0.10 | 4.17 | 17.46 | 19.26 | 19.27 | 5.26 
SUMT 5932.00 | 30.69 | 2.37 | 31.62 | 11.66 | 0.10 | 3.71 | 21.71 | 20.90 | 13.97 | 3.26 
LINRM | 6249.00 | 21.57 | 10.98 | 22.08 | 14.95 | 0.10 | 10.98 | 18.91 | 18.42 | 18.40 | 13.51 
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The comparison of the obtained results with [17] confirms the PSO algorithm used in this 
research for the truss optimization. 


4.2. Structural reliability analysis 


In this section, the reliability level of the truss system is determined to evaluate the difference 
between the available optimization approaches. The reliability analysis process is established 
based on the provided framework by Ghasemi and Nowak [18,19]. To do so, two types of the 
limit state functions based on the ultimate limit state function (ULSF). 


ULSF: gy = oy —= (30) 


where g,, and g, stand as the ultimate and service limit state functions, respectively. P is the 
element’s load; A represents the optimum value of the element’s cross-section. Tables 5 shows 
the reliability indices of each truss element, corresponding to the considered optimization 
approaches. It is worth mentioning that the reliability index is simply taken from FORM methods 
as shown in the following equation [20]. 


Ri Ans0j4 
Be = (31) 
[oR+03 


where R, and Q, are the nominal value of the structural resistance (nominal value of o,,for 
USLF and nominal value of aLfor SLSF) and nominal value of loads. Agand Agdenote the bias 
factor (mean over nominal value) of resistance and load. Also, ogand ogare the standard 
deviation of the resistance and load. The statistical parameter of both resistance and load can be 
taken from Ghasemi and Lee [21]. 


Table 5 
Reliability indices of each element corresponding to the given optimization approaches concerning ULSF. 


Method | Weight (Ib) 
OPTDYN 5472.00 
CONMIN 5563.00 


5620.08 
5613.84 


5719.00 
5725.00 
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As can be seen in Figure 7 and based on the load distribution analysis, elements Al, A4, A8 and 
A9 do not contributed to carry the applied loads. Also, since the truss has a two degree of 
indeterminists, we let two of the elements to exceed the allowable stress at a joint in such a way 
that elimination of those element cannot lead to the system instability. A5 to fail due to the high 
amount of applied tension. Therefore, in this illustrative example, elimination of two elements 
out of A5, A7, and A10 is permissible. Accordingly, GENETIC, Rajeev, GRP-UI, LINRM, and 
PSO are the only options which are resulting from a acceptable reliability indices greater than 
zero. Accordingly, a minimal reliability index of A2, A3, and A6 will be a dominant reliability 
index of a considered method. 


Table 6 
System reliability indices comparison. 
Method Weight (Ib) B 
GENETIC 5620.08 3.34 
Rajeev 5613.84 3.46 
GRP-UI 5727.00 1.53 
LINRM 6249.00 4.47 
PSO 5636.20 4.64 


As tabulated in Table 6. the highest reliability index is related to the PSO method. Although the 
second reliable approach is LINRM, the assigned cross-sections with correspondence to this 
method is the heaviest and the expensive option. 


5. Result and discussion 


In computational science, molecule swarm optimization (PSO) could be a computational method 
that optimizes an issue by iteratively attempting to make strides a candidate arrangement with 

respect to a given degree of quality. It understands an issue by having a populace of candidate 

arrangements, here named particles, and moving these particles around within the search-space 

concurring to basic numerical formulae over the particle's position and speed. Each particle's 
movement is affected by its local best-known position but is additionally guided toward the 
leading known positions within the search-space, which are overhauled as a way better positions 
are found by other particles. This is often anticipated to move the swarm toward the finest 

arrangements. 


PSO may be a metaheuristic because it makes few or no suspicions approximately the issue 
being optimized and can look very large spaces of candidate arrangements. Be that as it may, 
metaheuristics such as PSO don't ensure an ideal arrangement is ever found. Moreover, PSO 
does not utilize the slope of the issue being optimized, which implies PSO does not require that 
the optimization issue be differentiable as is required by classic optimization strategies such as 
slope plunge and quasi-newton strategies. 
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PSO has to been connected to multi-objective issues, in which the objective work comparison 
takes Pareto dominance into consideration when moving the PSO particles and non-dominated 
arrangements are put away to inexact the Pareto front. Multi-swarm optimization could be a 
variation of particle swarm optimization (PSO) based on the utilize of numerous sub-swarms 
rather than one (standard) swarm. The common approach in multi-swarm optimization is that 
each sub-swarm centers on a particular locale whereas a particular enhancement strategy chooses 
where and when to dispatch the sub-swarms. The multi-swarm system is particularly fitted for 
the optimization of multi-modal issues, where numerous (local) optima exist. 


6. Conclusions 


In this paper, a two-dimensional truss was analyzed with the finite element method and 
optimized using the PSO algorithm and then, the optimization results were compared with those 
of [17]. The analysis of comparisons suggests that the PSO algorithm leads to better results in 
terms of the higher convergence rate. Comparing the results with this reference shows confirms 
the PSO algorithm in this paper for truss optimization. It seems that the PSO algorithm (bird 
algorithm) is a good choice for the optimization of trusses and has an acceptably high rate of 
convergence. This method is dependent on the values of @, c:, and c2, and the values suggested in 
this paper can be used as the initial data. The number of particles (birds) is also important in this 
algorithm, but it is less important than the above factors and can be considered from 40 to 70. 


By the way, this problem opens two main question in structural optimizations. The first question 
is to how establish a framework to determine the reliability index of a system associated with the 
different degree of indeterminacies? And the second question is how to make a decision of 
optimization approach selection with consideration of both cost and reliability? The 
methodologies of addressing the responses to the above-mentioned questions are expected to be 
conducted in future studies. This, this paper attempts to shed more light to address the concern of 
the system reliability with various possible redundancy and target reliability selection associated 
with cost and reliability. 
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